GAPSO-Optimized Fuzzy PID Controller for Electric-Driven Seeding

To improve the seeding motor control performance of electric-driven seeding (EDS), a genetic particle swarm optimization (GAPSO)-optimized fuzzy PID control strategy for electric-driven seeding was designed. Since the parameters of the fuzzy controller were difficult to determine, two quantization factors were applied to the input of the fuzzy controller, and three scaling factors were introduced into the output of fuzzy controller. Genetic algorithm (GA) and particle swarm optimization (PSO) were combined into GAPSO by a genetic screening method. GAPSO was introduced to optimize the initial values of the two quantization factors, three scaling factors, and three characteristic functions before updating. The simulation results showed that the maximum overshoot of the GAPSO-based fuzzy PID controller system was 0.071%, settling time was 0.408 s, and steady-state error was 3.0693 × 10−5, which indicated the excellent control performance of the proposed strategy. Results of the field experiment showed that the EDS had better performance than the ground wheel chain sprocket seeding (GCSS). With a seeder operating speed of 6km/h, the average qualified index (Iq) was 95.83%, the average multiple index (Imult) was 1.11%, the average missing index (Imiss) was 3.23%, and the average precision index (Ip) was 14.64%. The research results provide a reference for the parameter tuning mode of the fuzzy PID controller for EDS.


Introduction
The quality of seeding greatly depends on the rotating speed of the metering disc driven by the seeding motor. Motor control strategies mainly include PID control, adaptive control, vector control, and intelligent control. The intelligent control has a blocked hierarchical control organization structure, which is convenient for handling a large amount of information, stored knowledge ability, and reasoning ability. Intelligent control carried out real-time control of the dynamic environment of the seeding motor movement process according to the characteristics of nonlinear, time-varying, dynamic, and complex uncertain control objects under the steady state, when driving the seeding motor speed regulation in the electric-driven seeding (EDS) [1][2][3]. Intelligent control strategies mainly include: fuzzy control [4][5][6][7][8], neural network [9][10][11][12], and the heuristic optimization algorithm [13][14][15][16][17]. Among them, the heuristic optimization algorithm has a significant role in tuning the parameters of the controller. The stronger the global optimization ability of the algorithm, the better the robustness of the controller adjusted by the algorithm. In recent years, scholars in various fields have also improved and applied the heuristic optimization algorithm and achieved remarkable achievements [18][19][20].
Zhai et al. [21] designed a control system for precision seed-metering device. The system adjusts the relationship between the operating speed of the seeder and the seeding motor speed, according to seed spacing and seeding rate of the seeder. Zhang et al. [22] designed an electronic control seeding system for maize. The desired rotating speed and actual rotating speed of the seed-metering device were taken as inputs, and the PID

Modeling for Seeding Motor
The research object was the air suction mung bean-metering device. Assuming that there was m seed cells in the metering disc of the seed-metering device, the time difference between two adjacent mung bean seeds falling is as follows: where ∆t is the falling time difference between two mung bean seeds, s; n is the rotating speed of the seed-metering device, r/min. The equation of mung bean seed spacing is Z = 277.78v∆t = 1.67 × 10 4 v nm (2) where Z is mung bean seed spacing, mm; v is the operating speed of seeder, km/h. The seed cells number m is 27 in the mung bean seed-metering disc in the selected seed-metering device, so the following equation is obtained from Equation (2).
According to the technical specification for the Heilongjiang Province local standard Technical Regulations for Mechanized Production of Mung Bean (DB23/T 2767-2020), when mung beans are sown in a single row by a tractor with more than 25 horsepower, the seed spacing is 80 mm~100 mm. Take 100 mm in study and substitute the plant spacing value into Equation (3) to obtain: n = 6.19v (4) In order to establish the simulation system for designed intelligent control strategy of the seeding motor, it was necessary to establish a mathematical model of the seeding motor. JGB37-3650-1280 DC brushless reducer motor produced by ASLONG Company was selected as the seeding motor for research. The seeding motor parameters are shown in Table 1. The mechanical characteristics equation of brushless DC motor is the same as that of separately excited DC motor. Because the fluctuation of back-EMF and torque of DC brushless motor is relatively large, the concepts of average torque and average back-EMF are used in the mechanical characteristics equation. According to the principle of brushless DC motor, the voltage-balance equation of brushless DC motor is obtained: where U 0 is armature voltage, V; L is armature inductance, H; I is armature current, A; R is armature resistance, Ω; C e is the back-EMF coefficient, V·(rad·s −1 ) −1 ; n is the seeding motor speed, r/min. The torque-balance equation of brushless DC motor is as follows: where T e is the electromagnetic torque, g·cm; J is the moment of inertia, kg·m 2 ; T f is the load torque, g·cm. According to Equations (5) and (6), the differential equation of DC brushless motor can be calculated as follows: where T d is the electromagnetic time constant, ms; T m is the mechanical time constant, ms. By Laplace transform from Equation (7), the transfer function of brushless DC motor can be obtained as follows: The T d and T m are calculated as follows: where K t is the torque constant, N·m·A −1 .
The T d and T m can be calculated by Equations (9) and (10). The seeding motor parameters T m is 0.017, T d is 0.855, and C e is 0.192. They are substituted into Equation (8) to obtain the seeding motor mathematical model:

Establishment of Adaptive Fuzzy Controller
Compared with the general fuzzy system, the adjustment process of the error and error rate of the whole system was controlled by the same fuzzy reasoning rules. The change range of error e and error rate e c of the system was defined as the domain on the fuzzy set. According to the set fuzzy inference rules, the system can adaptively modify k p , k i and k d . The modified parameters could be obtained by deblurring the obtained results ∆k p , ∆k i , and ∆k d , which were the three outputs of the fuzzy controller. After substituting them into Equations (12)- (14), respectively, the three modified characteristic functions k p , k i , and k d can be obtained by calculation.
where k p(k−1) , k i(k−1) , and k d(k−1) are the three PID characteristic functions before updating; ∆k p(k) , ∆k i(k) , and ∆k d(k) are the correction parameters of PID characteristic functions; k p(k) , k i(k) , and k d(k) represent the updated PID characteristic functions. Then, k p(k) , k i(k) , and k d(k) were placed into the transfer function of PID and the online automatic correction of PID parameters was completed. In order to establish a fuzzy PID controller, two input and three output structures were chosen for the fuzzy controller. Mamdani method was selected as the fuzzy inference method, and the membership function of the fuzzy controller is shown in Figure 1.  According to the membership assignment rules of each fuzzy subset and the fuzzy control model of each parameter, the PID parameters could be modified online by adopting fuzzy control rules. The fuzzy matrix table of PID modification parameters Δkp, Δki, and Δkd was designed by using fuzzy synthetic reasoning. The fuzzy inference rules were 49 kinds of conventional fuzzy inference rules. According to the membership degree assignment table of each fuzzy subset and the fuzzy control model of each parameter, the PID parameters can be modified online by using the fuzzy control rules. The fuzzy matrix table of PID controller correction parameters is designed by using fuzzy synthesis reasoning, and the designed PID correction parameters are studied Δkp, Δki, and Δkd. The 49 fuzzy inference rules are shown in Table 2. The fuzzy surfaces, corresponding to the modified parameters, are shown in Figure 2. According to the membership assignment rules of each fuzzy subset and the fuzzy control model of each parameter, the PID parameters could be modified online by adopting fuzzy control rules. The fuzzy matrix table of PID modification parameters ∆k p , ∆k i , and ∆k d was designed by using fuzzy synthetic reasoning. The fuzzy inference rules were 49 kinds of conventional fuzzy inference rules. According to the membership degree assignment table of each fuzzy subset and the fuzzy control model of each parameter, the PID parameters can be modified online by using the fuzzy control rules. The fuzzy matrix table of PID controller correction parameters is designed by using fuzzy synthesis reasoning, and the designed PID correction parameters are studied ∆k p , ∆k i , and ∆k d . The 49 fuzzy inference rules are shown in Table 2. The fuzzy surfaces, corresponding to the modified parameters, are shown in Figure 2. control model of each parameter, the PID parameters could be modified online by adopting fuzzy control rules. The fuzzy matrix table of PID modification parameters Δkp, Δki, and Δkd was designed by using fuzzy synthetic reasoning. The fuzzy inference rules were 49 kinds of conventional fuzzy inference rules. According to the membership degree assignment table of each fuzzy subset and the fuzzy control model of each parameter, the PID parameters can be modified online by using the fuzzy control rules. The fuzzy matrix table of PID controller correction parameters is designed by using fuzzy synthesis reasoning, and the designed PID correction parameters are studied Δkp, Δki, and Δkd. The 49 fuzzy inference rules are shown in Table 2. The fuzzy surfaces, corresponding to the modified parameters, are shown in Figure 2.   In the Simulink module, the quantization coefficients k e and k ec can be introduced to the inputs e and e c as the quantization factors of the fuzzy controller. The definition of the quantization factors k e and k ec is: n e e max , k ec = n ec ec max (15) where k e is the quantization factor of error, k ec is the quantization factor of error change rate, n e is the fuzzy series of error, n ec is the fuzzy series of error change rate, e max is the maximum error value, and ec max is the error value of error change rate. PID parameter calibration is carried out by formulating fuzzy rules of the fuzzy controller in advance. According to Equations (12)- (14), the output ∆k p , ∆k i , and ∆k d adaptively adjusts three parameters k p , k i , and k d in the PID controller; in the Simulink module, ∆k p , ∆k i , and ∆k d introduced the scaling factors dk p , dk i , and dk d as the scaling factors of the fuzzy controller. The definition of the scaling factors dk p , dk i , and dk d is: where dk p is the scaling factor of the proportional correction parameter; dk i is the scaling factor of the differential correction parameter; dk d is the scaling factor of the integral correction parameter; k pmax , k imax , and k dmax are the maximum change values of the three characteristic parameters of the PID controller; and N is the fuzzy output series of the fuzzy controller. It can be found clearly from Figures 1 and 2, the fuzzy domain of fuzzy input and fuzzy output was taken as {−3, −2, −1, 0, 1, 2, 3}. Because there was no fixed rules for the selection of the fuzzy domain and the setting process of two inputs and three outputs, the fuzzy domain is cumbersome. The essence of fuzzy control was to change the fuzzy series under the preset fuzzy rules through fuzzy input and fuzzy output. The quantitationscaling coefficient could be introduced into the fuzzy input and fuzzy output, respectively, and the fuzzy series could be adjusted by optimizing the coefficient, so as to achieve the best control effect. Figure 3 shows the simulink model of fuzzy PID control system. The control system takes the theoretical speed of the seed metering device as the step input and makes the difference between the input and the feedback through the comparator to obtain the system error e. The system error is differentiated and the system error change rate ec is obtained. e and ec are the inputs of the fuzzy controller, and the three correction values of the PID controller can be obtained through the fuzzy rules. The correction value and the set value can be added to obtain the three modified characteristic functions, the corrected characteristic function is multiplied by the error, the integral characteristic function is multiplied by the error integral value, and differential characteristic function is multiplied by error differential values. The above three values are characteristic functions that affect the control effect of the PID controller, and then the PID controller outputs the control amount to the seeding motor until the target value is reached. Therefore, the introduction of scaling factors and quantization factors can directly change the correction ability of the characteristic functions of the PID controller. Quantization coefficients k e and k ec were used as quantization factors of inputs e and e c of the fuzzy controller. To the fuzzy outputs ∆k p , ∆k i , and ∆k d , the scaling coefficients dk p , dk i , and dk d were used as the scaling factor of the fuzzy controller. According to Z-N method, the values of three characteristic functions before the updating of the PID controller were k p = 3.5, k i = 0.9, and k d = 0.08. All quantization factors and scaling factors were set to 1. After experiment, the results of the EDS system were as follows: the maximum overshoot was 6.399%, settling time was 3.31 s, and steady-state error was −0.5287.  Fuzzy PID had obvious disadvantages, that there was no general setting rule for membership functions. The fuzzy rules and membership functions after setting the fuzzy series of the fuzzy controller could not be adjusted. The process of setting the domain and fuzzy rules was complicated; the process of resetting the domain was complex; the ideal control effect might not be obtained; and maximum overshoot, settling time, and steadystate error are difficult to achieve optimum.
In the seeding, because of the uneven terrain, the operating speed will be often disturbed, so the settling time of the system has a great influence on the results. In the seeding, with the increase in operating time, the steady-state error accumulation of the system will become larger and larger, which is the total seeding rate deviation caused by the total stroke deviation of the seeding motor, so the system should also keep small steady-state error.

GAPSO-Optimized Fuzzy PID Controller for EDS
PSO is a swarm intelligence algorithm based on birds flocking behavior. Let NP be the number of particles in PSO. A particle Pi has position Xi and velocity Vi. The fitness function is used to evaluate each particle for checking the quality of the solution. Initially, each particle is assigned with a random position and velocity values. Each particle computes its own best called Pbesti and global best called Gbest for every iteration [30]. To reach the global best solution, it uses its personal and global best to update the velocity Vi Fuzzy PID had obvious disadvantages, that there was no general setting rule for membership functions. The fuzzy rules and membership functions after setting the fuzzy series of the fuzzy controller could not be adjusted. The process of setting the domain and fuzzy rules was complicated; the process of resetting the domain was complex; the ideal control effect might not be obtained; and maximum overshoot, settling time, and steady-state error are difficult to achieve optimum.
In the seeding, because of the uneven terrain, the operating speed will be often disturbed, so the settling time of the system has a great influence on the results. In the seeding, with the increase in operating time, the steady-state error accumulation of the system will become larger and larger, which is the total seeding rate deviation caused by the total stroke deviation of the seeding motor, so the system should also keep small steady-state error.

GAPSO-Optimized Fuzzy PID Controller for EDS
PSO is a swarm intelligence algorithm based on birds flocking behavior. Let N P be the number of particles in PSO. A particle P i has position X i and velocity V i . The fitness function is used to evaluate each particle for checking the quality of the solution. Initially, each particle is assigned with a random position and velocity values. Each particle computes its own best called Pbest i and global best called Gbest for every iteration [30]. To reach the global best solution, it uses its personal and global best to update the velocity V i and position X i using the following Equations (17) and (18).
where w is the inertia weight, c 1 , c 2 are acceleration coefficients and χ 1 , χ 2 are randomly generated values. After getting a new updated position, the particle evaluates the fitness function and updates Pbest i as well as Gbest from Equations (19) and (20).
The GA is based on the principles of natural selection. It mimics the process of survival of the fittest principle in nature by trying to maximize the fitness function. The population, which represents the optimization variable sets, is updated after each learning cycle through three evolutionary processes, i.e., selection, crossover, and mutation [31].
PSO has the memory ability, and the optimal filial-population particles can be preserved, and per-generation memory will not be forgotten because of population evolution. GA does not have memory ability, and the per-generation memory will be forgotten because of population evolution. In terms of convergence, GA has a perfect convergence analysis method, and estimates the convergence rate. However, the convergence analysis on PSO is still weakness. Although there have been a simplified version of convergence analysis, the transformation from certainty to randomness needs to be further studied. In addition, unlike GA, PSO has no evolution tools, such as crossover or mutation [31]. Both GA and PSO have disadvantages, as genetic screening combination of them can not only enhance the memory ability but also enhance the convergence ability.
The study adopts the genetic screening mechanism and proposes a GAPSO combined GA and PSO. Compared with the other two algorithms, its memory ability is enhanced with the optimal solution preservation of PSO. At the same time, it can also increase the global optimization probability. The combination steps of GA and PSO are as follows: Step 1. Initialize the particle swarm and set the initial parameters. Such as population size, crossover probability, mutation probability, particle dimension, constraint matrix.
Step 2. Calculate the fitness value of the particle swarm using the objective function, determine the optimal position of the individual particle and particle swarm, update the optimal particle, and calculate the optimal position of individual particles (Pbest i ) of the particle swarm.
Step 3. Judge whether the ultimate objective of the optimization process is achieved. The ultimate objective is usually preset as the maximum number of iterations and the fitness value of the GAPSO reaching the standard. If the ultimate objective is achieved, the iterative process will be stopped, otherwise, Step 4 will be continued.
Step 4. Update particle velocity, particle position, and optimal value. The roulette wheel selection was used to select the particle swarm. It was decided to keep three quarters Sensors 2022, 22, 6678 8 of 18 of the particle swarm and eliminate a quarter of the particle swarm. The probability distribution equation of the selection probability P i is as follows: where P i is the selection probability of particle i; n is the population size; f(x i ) is the single particle fitness of particle x i ; f(x i ) is the fitness of all particles.
Step 5. Crossover in the intermediate generation particles. The new generation of individuals is generated by adaptive arithmetic uniform crossover with the crossover probability P c of 80%. The calculation equation of P c is as follows: where P c is the crossover probability; f max is the maximum fitness of particle swarm; f avg is the average fitness of particle swarm; f is the maximum fitness of two crossover particles; k 1 and k 2 are constant.
Step 6. The non-consistent mutation is used to mutate part of all particles according to the mutation probability P m . Because it needs to conform to the natural law, P m is generally set to be less than 0.05. In order to avoid falling into local optimum. P m is 0.1 in this study. The equation for calculating P m is as follows: where P m is the mutation probability; k 3 and k 4 are constant.
Step 7. Update the particle swarm. The Pbest i is updated by calculating, and then return to Step 4. When the preset maximum iteration is reached, Gbest is output, and the GAPSO stops.
The process of setting control parameters by GAPSO is shown in Figure 4. The input signal was the operating speed, and the desired rotating speed of seedmetering device was calculated from Equation (3). The seeding motor driven the seedmetering device to rotate, and the feedback value of the system was the rotor speed value of the seeding motor. The difference between the desired rotating speed of the seed-metering device and the actual rotating speed value was calculated as the error e, and the error change rate e c was the differential of the error e. The flowchart of the GAPSO-optimized fuzzy PID controller for EDS is shown in Figure 5.
The integral of time multiplied by the absolute value of error criterion (ITAE) is taken as an evaluation index in GAPSO. The optimized parameters are the initial values of three characteristic functions before updating in the PID controller, two quantization factors and three scaling factors in the fuzzy controller. Thus, the step response curve of the optimality is obtained in set iteration number.
The system designed according to the index has less oscillation and selectivity to parameters. The ITAE was used in the objective function in the study, which synthesizes the total running time and the absolute value of the system steady-state error; the smaller its value, the better the performance of the system, and the equation is: where t is the system running time; e(t) is the system steady-state error. ITAE index could evaluate the static and dynamic performance of the system, such as rapidity, accuracy, and stability. This index could achieve remarkable optimization effect on the steady-state error and settling time of the seeding motor. Step 7. Update the particle swarm. The Pbesti is updated by calculating, and then return to Step 4. When the preset maximum iteration is reached, Gbest is output, and the GAPSO stops.
The process of setting control parameters by GAPSO is shown in Figure 4. The input signal was the operating speed, and the desired rotating speed of seedmetering device was calculated from Equation (3). The seeding motor driven the seedmetering device to rotate, and the feedback value of the system was the rotor speed value of the seeding motor. The difference between the desired rotating speed of the seed-metering device and the actual rotating speed value was calculated as the error e, and the error change rate ec was the differential of the error e. The flowchart of the GAPSO-optimized fuzzy PID controller for EDS is shown in Figure 5. The integral of time multiplied by the absolute value of error criterion (ITAE) is taken as an evaluation index in GAPSO. The optimized parameters are the initial values of three characteristic functions before updating in the PID controller, two quantization factors and three scaling factors in the fuzzy controller. Thus, the step response curve of the optimality is obtained in set iteration number.
The system designed according to the index has less oscillation and selectivity to parameters. The ITAE was used in the objective function in the study, which synthesizes the

Simulation Experiences and Results
Because there are eight parameters to be optimized, the overlong simulation running time should be avoided. Before the simulation, the iteration range of these eight parameters is set according to the previous debugging experience of the fuzzy PID controller. The iteration ranges of the three characteristic functions before updating in the PID controller are: In the coding process, the fuzzy PID controller adopts the form of two inputs and three outputs, and each input and output corresponds to seven fuzzy subsets. The fuzzy subsets from negative to positive are represented by integers from 1 to 7. The fuzzy control rule table is transformed into a two-dimensional matrix. The two-dimensional matrix is straightened and arranged to form chromosome individuals, and the length of chromosome is 7 2 × 3 = 147.
In the GA, the crossover probability P c is set to 0.8 (80%), the mutation probability P m is set to 0.1 (10%), the chromosomes number is 1, the chromosomes length is 147, the population size n is set to 20, and the maximum number of iterations is 100.
In the PSO, the inertia coefficient ω is set to 1, the acceleration coefficients C 1 and C 2 are 1.4945, the population number is 10, the population size n is 20, and the maximum number of iterations is 100.
In the GAPSO, the crossover probability P c is set to 0.8 (80%), the mutation probability P m is set to 0.1 (10%), the chromosomes number is 1, the chromosome length is 147, the population size n is set to 20, the inertia coefficient ω is set to 1, the acceleration coefficients C 1 and C 2 are set to 1.4945, and the maximum number of iterations is 100.
The Matlab2020b is used as the simulation environment of experiment. The performances of GA, PSO, and GAPSO on optimizing fuzzy PID controller parameters were compared, and the advantages and disadvantages of the designed control strategy were discussed. The optimization performance data table of the three control strategies are shown in Table 3. The unit step response curves optimized by the three control strategies are shown in Figure 6. The fitness function control curves are shown in Figure 7. In the GA, the crossover probability Pc is set to 0.8 (80%), the mutation probability Pm is set to 0.1 (10%), the chromosomes number is 1, the chromosomes length is 147, the population size n is set to 20, and the maximum number of iterations is 100.
In the PSO, the inertia coefficient ω is set to 1, the acceleration coefficients C1 and C2 are 1.4945, the population number is 10, the population size n is 20, and the maximum number of iterations is 100.
In the GAPSO, the crossover probability Pc is set to 0.8 (80%), the mutation probability Pm is set to 0.1 (10%), the chromosomes number is 1, the chromosome length is 147, the population size n is set to 20, the inertia coefficient ω is set to 1, the acceleration coefficients C1 and C2 are set to 1.4945, and the maximum number of iterations is 100.
The Matlab2020b is used as the simulation environment of experiment. The performances of GA, PSO, and GAPSO on optimizing fuzzy PID controller parameters were compared, and the advantages and disadvantages of the designed control strategy were discussed. The optimization performance data table of the three control strategies are shown in Table 3. The unit step response curves optimized by the three control strategies are shown in Figure 6. The fitness function control curves are shown in Figure 7.     According to Figure 6 and Table 3, compared with fuzzy PID control, the various performance indexes have been improved by three intelligent control strategies. Compared with the other two intelligent control strategies, GAPSO-optimized fuzzy PID controller had advantages in the maximum overshoot, settling time, and steady-state error. The maximum overshoot was 0.071%, the settling time was 0.408 s, the steady-state error was 3.0693e −5 , and that was no oscillation after reaching steady-state. Although the control strategy of optimizing fuzzy PID parameters by GA and PSO obviously improved the performance of the system, compared with fuzzy PID, they was still a poor performance in optimizing maximum overshoot and settling time, and there was still a small amount of oscillation after reaching steady-state.
According to Figure 7 and Table 3, the GA-optimized fuzzy PID controller reached the optimum in the 41st iteration, and the fitness value was 0.7531. The PSO-optimized fuzzy PID controller reached the optimum in the 48th iteration, and the fitness value was 0.6311. The GAPSO-optimized fuzzy PID controller reached the optimum in the 65th iteration, and the fitness value was 0.5575. Compared with the other two intelligent control strategies, GAPSO needed more iterations to optimize the fuzzy PID controller. However, compared with the other two intelligent algorithms, GAPSO had the smaller fitness value optimized and the better system performance.
With the GAPSO-optimized fuzzy PID controller, at the 65th iteration, the value of ke was 6.5752. At the 65th iteration, the value of kec was 7.5865. At the 5th iteration, the value of dkp was 0.0060. At the 5th iteration, the value of dki was 14.5898. At the 4th iteration, the value of dkd was 0.0020. At the 5th iteration, the value of kp was 30.1425. At the 65th iteration, the value of ki was 14.2584. At the 40th iteration, the value of kd was 0.3333.
In the optimization process of quantization factors and three characteristic functions before updating, the GA-optimized fuzzy PID controller had more advantages than the other two control strategies and could faster iterate the optimal value. In the optimization process of scaling factors, the PSO-optimized fuzzy PID controller had more advantages than the other two control strategies and could reach the optimal value faster. Although the GAPSO-optimized fuzzy PID controller has many iterations, it had smaller fitness value and better control performance. Therefore, the intelligent control strategy was selected as the EDS control strategy in this study. According to Figure 6 and Table 3, compared with fuzzy PID control, the various performance indexes have been improved by three intelligent control strategies. Compared with the other two intelligent control strategies, GAPSO-optimized fuzzy PID controller had advantages in the maximum overshoot, settling time, and steady-state error. The maximum overshoot was 0.071%, the settling time was 0.408 s, the steady-state error was 3.0693e −5 , and that was no oscillation after reaching steady-state. Although the control strategy of optimizing fuzzy PID parameters by GA and PSO obviously improved the performance of the system, compared with fuzzy PID, they was still a poor performance in optimizing maximum overshoot and settling time, and there was still a small amount of oscillation after reaching steady-state.
According to Figure 7 and Table 3, the GA-optimized fuzzy PID controller reached the optimum in the 41st iteration, and the fitness value was 0.7531. The PSO-optimized fuzzy PID controller reached the optimum in the 48th iteration, and the fitness value was 0.6311. The GAPSO-optimized fuzzy PID controller reached the optimum in the 65th iteration, and the fitness value was 0.5575. Compared with the other two intelligent control strategies, GAPSO needed more iterations to optimize the fuzzy PID controller. However, compared with the other two intelligent algorithms, GAPSO had the smaller fitness value optimized and the better system performance.
With the GAPSO-optimized fuzzy PID controller, at the 65th iteration, the value of k e was 6.5752. At the 65th iteration, the value of k ec was 7.5865. At the 5th iteration, the value of dk p was 0.0060. At the 5th iteration, the value of dk i was 14.5898. At the 4th iteration, the value of dk d was 0.0020. At the 5th iteration, the value of k p was 30.1425. At the 65th iteration, the value of k i was 14.2584. At the 40th iteration, the value of k d was 0.3333.
In the optimization process of quantization factors and three characteristic functions before updating, the GA-optimized fuzzy PID controller had more advantages than the other two control strategies and could faster iterate the optimal value. In the optimization process of scaling factors, the PSO-optimized fuzzy PID controller had more advantages than the other two control strategies and could reach the optimal value faster. Although the GAPSO-optimized fuzzy PID controller has many iterations, it had smaller fitness value and better control performance. Therefore, the intelligent control strategy was selected as the EDS control strategy in this study.

Speed Accuracy Experiences of Seed Metering Device
In order to explore the speed accuracy of the three intelligent electric-driven mung bean air suction precision seed metering device control systems studied; the actual speed of the seed metering device at the theoretical seed spacing of 100 mm; the operation speed of 3 km/h, 5 km/h, 7 km/h, and 9 km/h; the speed of 30 s after the speed of the seed metering motor is stabilized; and the average speed N of the seed metering device is calculated through the DT-2234A digital photo tachometer (accuracy: ±0.1 rpm) developed by Lutron company. The speed of each group is repeated five times (N 1~N5 ), and the variation coefficient (CV), mean, standard deviation (SD), and confidence interval (CI, 99% confidence) of the actual speed of the seed metering device is calculated. In order to compare the actual control effect of three intelligent control strategies. Figure 8 is the seed metering device speed accuracy experiment.

Speed Accuracy Experiences of Seed Metering Device
In order to explore the speed accuracy of the three intelligent electric-driven mung bean air suction precision seed metering device control systems studied; the actual speed of the seed metering device at the theoretical seed spacing of 100 mm; the operation speed of 3 km/h, 5 km/h, 7 km/h, and 9 km/h; the speed of 30 s after the speed of the seed metering motor is stabilized; and the average speed N of the seed metering device is calculated through the DT-2234A digital photo tachometer (accuracy: ±0.1 rpm) developed by Lutron company. The speed of each group is repeated five times (N1~N5), and the variation coefficient (CV), mean, standard deviation (SD), and confidence interval (CI, 99% confidence) of the actual speed of the seed metering device is calculated. In order to compare the actual control effect of three intelligent control strategies. Figure 8 is the seed metering device speed accuracy experiment. Table 4 is the verification experiment table of the speed accuracy of the seed metering device.   It can be seen from Table 4 that the mung bean EDS control system based on the GAPSO-Fuzzy-PID algorithm has a good precision effect of the seed metering device speed, good uniformity of the speed, and the variation coefficient of the seed metering device speed is less than 6.00%. When the operation speed is 9 km/h, the maximum variation coefficient of the seed metering device actual speed reaches 5.90%; when the operating speed is 7 km/h, the minimum variation coefficient of the seed metering device actual speed is 2.99%. According to the standard deviation, it can be shown that GAPSO fuzzy PID has higher speed accuracy of the seed metering device, compared with the other two intelligent control systems, especially in the high-speed seeding state. At 7 km/h, the average speed of GAPSO fuzzy PID is the closest to the theoretical value, and the difference is only 0.02. In addition, the deviation rate of the seed metering device actual speed is small, and the speed deviation rate is less than 12%. When the operation speed is 9 km/h, N 4 reaches the maximum deviation speed, and the speed deviation rate is 11.31%; when the operation speed is 7 km/h, N 5 reaches the minimum deviation speed, and the speed deviation rate is 0.69%. Therefore, GAPSO-Fuzzy-PID with higher rotational speed accuracy is selected for field experiments.

Field Experiments and Results
The field experiment was conducted on April 20th, 2021 in the Heilongjiang Province Green Grassland Pasture experimental field, zhonglv-5 mung bean hybrid was selected as experiment materials, and two-rows air-suction precision seeder reformed by Heilongjiang Bayi Agricultural University was used for the experiment. The left seeding unit was the ground wheel chain sprocket seeding (GCSS), and the right seeding unit was the EDS. The field experiment is shown in Figure 9. It can be seen from Table 4 that the mung bean EDS control system based on the GAPSO-Fuzzy-PID algorithm has a good precision effect of the seed metering device speed, good uniformity of the speed, and the variation coefficient of the seed metering device speed is less than 6.00%. When the operation speed is 9 km/h, the maximum variation coefficient of the seed metering device actual speed reaches 5.90%; when the operating speed is 7 km/h, the minimum variation coefficient of the seed metering device actual speed is 2.99%. According to the standard deviation, it can be shown that GAPSO fuzzy PID has higher speed accuracy of the seed metering device, compared with the other two intelligent control systems, especially in the high-speed seeding state. At 7 km/h, the average speed of GAPSO fuzzy PID is the closest to the theoretical value, and the difference is only 0.02. In addition, the deviation rate of the seed metering device actual speed is small, and the speed deviation rate is less than 12%. When the operation speed is 9 km/h, N4 reaches the maximum deviation speed, and the speed deviation rate is 11.31%; when the operation speed is 7 km/h, N5 reaches the minimum deviation speed, and the speed deviation rate is 0.69%. Therefore, GAPSO-Fuzzy-PID with higher rotational speed accuracy is selected for field experiments.

Field Experiments and Results
The field experiment was conducted on April 20th, 2021 in the Heilongjiang Province Green Grassland Pasture experimental field, zhonglv-5 mung bean hybrid was selected as experiment materials, and two-rows air-suction precision seeder reformed by Heilongjiang Bayi Agricultural University was used for the experiment. The left seeding unit was the ground wheel chain sprocket seeding (GCSS), and the right seeding unit was the EDS. The field experiment is shown in Figure 9. The length of the plot in the experiment area was 60 m, and the preparation areas at both ends were 20 m, respectively, 120 groups of stable continuous operation data were collected in the middle section 20 m, repeated three times, and the average value of the results were calculated. Fixed parameters of the experiment: the wind pressure of the blower was 2.5 kpa, and the starting angle of seeding was 20°. According to the Chinese national standard test method of single seed (Precision) seeder (GB/T 6973-2005). The qualified index (Iq), multiple index (Imult), missing index (Imiss), and precision index (Ip) were selected as the index of seeding performance. According to the method, when the theoretical seeding spacing is 100 mm, the qualified range of seeding spacing is 50-150 mm. Under different operating speeds of 3 km/h, 6 km/h, 9 km/h, and 12 km/h, the Iq, Imult, Imiss, and Ip between GCSS and EDS were analyzed and compared, so as to verify the superiority of the control strategy studied. The experiment statistical results are shown in Table 5. The length of the plot in the experiment area was 60 m, and the preparation areas at both ends were 20 m, respectively, 120 groups of stable continuous operation data were collected in the middle section 20 m, repeated three times, and the average value of the results were calculated. Fixed parameters of the experiment: the wind pressure of the blower was 2.5 kpa, and the starting angle of seeding was 20 • . According to the Chinese national standard test method of single seed (Precision) seeder (GB/T 6973-2005). The qualified index (I q ), multiple index (I mult ), missing index (I miss ), and precision index (I p ) were selected as the index of seeding performance. According to the method, when the theoretical seeding spacing is 100 mm, the qualified range of seeding spacing is 50-150 mm. Under different operating speeds of 3 km/h, 6 km/h, 9 km/h, and 12 km/h, the I q , I mult , I miss , and I p between GCSS and EDS were analyzed and compared, so as to verify the superiority of the control strategy studied. The experiment statistical results are shown in Table 5. According to Table 5, various indexes had been promoted and improved by EDS in the field experiment, compared with the GCSS. According to the technical conditions of single seed (precision) seeder (JB/T 10293-2013), which was seed spacing ≤ 100 mm under, the I q ≥ 60%, the I mult ≤ 30%, the I miss ≤ 15%, and the I p ≤ 40%.
According to Figure 10, In each operation speed, each seeding performance of EDS was better than that of GCSS, when the operating speed was 6 km/h, both the GCSS and the EDS achieved the best seeding performance, but EDS compared with the GCSS, the I q of EDS was increased by 9.16%, the I mult was increased by 1.11%, the I miss was increased by 7.88%, and the I p was increased by 8.44%. At the operating speed of 12 km/h, the I miss of three GCSS experiments were higher than 15%. It showed that under the high-speed seeding operation, the GCSS had been difficult to meet the Chinese national standard, but the I miss of EDS could meet the Chinese national standard.  According to Table 5, various indexes had been promoted and improved by EDS in the field experiment, compared with the GCSS. According to the technical conditions of single seed (precision) seeder (JB/T 10293-2013), which was seed spacing ≤ 100 mm under, the Iq ≥ 60%, the Imult ≤ 30%, the Imiss ≤ 15%, and the Ip ≤ 40%.
According to Figure 10, In each operation speed, each seeding performance of EDS was better than that of GCSS, when the operating speed was 6 km/h, both the GCSS and the EDS achieved the best seeding performance, but EDS compared with the GCSS, the Iq of EDS was increased by 9.16%, the Imult was increased by 1.11%, the Imiss was increased by 7.88%, and the Ip was increased by 8.44%. At the operating speed of 12 km/h, the Imiss of three GCSS experiments were higher than 15%. It showed that under the high-speed seeding operation, the GCSS had been difficult to meet the Chinese national standard, but the Imiss of EDS could meet the Chinese national standard.  The electric-driven mung bean air suction precision seed metering device control system based on GAPSO-Fuzzy-PID can effectively reduce the missing index and precision index of the original seed metering device control system, reduce the system's false rate, make the seed spacing more uniform, and meet the qualified range of each index in the national standard. However, there was no significant increase in the multiple index. The distance distribution of GCSS under each index is shown in Figure 11. The distance distribution of intelligent EDS under each index is shown in Figure 12. The electric-driven mung bean air suction precision seed metering device control sys tem based on GAPSO-Fuzzy-PID can effectively reduce the missing index and precision index of the original seed metering device control system, reduce the system's false rate make the seed spacing more uniform, and meet the qualified range of each index in the national standard. However, there was no significant increase in the multiple index. The distance distribution of GCSS under each index is shown in Figure 11. The distance dis tribution of intelligent EDS under each index is shown in Figure 12. It can be seen from Figures 11 and 12 that, compared with the GCSS and the intelli gent EDS, the intelligent EDS can significantly reduce the missing index and the coefficien of variation in the mung bean seeding, which can make the dispersion of the actual seeded mung beans closer to the theoretical seed spacing line. This ensures that the mung bean seeds between the ridges are evenly aligned and the number is consistent, which is con venient for the follow-up field management of farmers. It can be seen from Figures 11 and 12 that, compared with the GCSS and the intelligent EDS, the intelligent EDS can significantly reduce the missing index and the coefficient of variation in the mung bean seeding, which can make the dispersion of the actual seeded mung beans closer to the theoretical seed spacing line. This ensures that the mung bean seeds between the ridges are evenly aligned and the number is consistent, which is convenient for the follow-up field management of farmers.

Discussions
(1) In the process of parameter optimization of the fuzzy PID controller, GAPSO had more iterations and longer off-line optimization time. Therefore, we should consider designing a new combination method of GA and PSO. The intelligent optimization algorithm with convergence and memory function could be used to shorten the off-line optimization time and improve the optimization efficiency.
(2) With the operating speed of 12 km/h, the I miss of the EDS was obviously higher than that in medium-and low-speed operation. It did not obviously change the drawback of the increase sharply of the I miss in the high-speed seeding operation of the GCSS. After removing a number of unstable factors, such as chain-driven, this might be related to the designed defects of the air-suction seed-metering device, and it might not be able to achieve high-speed seeding under GCSS.

Discussions
(1) In the process of parameter optimization of the fuzzy PID controller, GAPSO had more iterations and longer off-line optimization time. Therefore, we should consider de signing a new combination method of GA and PSO. The intelligent optimization algo rithm with convergence and memory function could be used to shorten the off-line opti mization time and improve the optimization efficiency.
(2) With the operating speed of 12 km/h, the Imiss of the EDS was obviously higher than that in medium-and low-speed operation. It did not obviously change the drawback of the increase sharply of the Imiss in the high-speed seeding operation of the GCSS. After removing a number of unstable factors, such as chain-driven, this might be related to the designed defects of the air-suction seed-metering device, and it might not be able to achieve high-speed seeding under GCSS.

Conclusions
This paper has presented a GAPSO-optimized Fuzzy PID controller for electric driven seeding. Through simulation experiments, speed accuracy experiments of seed me tering device, and field experiments, it is verified that the proposed GAPSO is better than the traditional GA and PSO in the parameter tuning of the fuzzy controller; compared with the other two algorithms, the GAPSO fused by genetic screening mechanism has stronger global optimization ability under the same iteration time.; Macroscopically, the rotational speed accuracy of the seeding motor driven by the three different intelligen controllers of GAPSO-Fuzzy-PID is also higher, and the field experiments results are ex cellent, which provides the basic conditions for the high-density planting of crops.

Conclusions
This paper has presented a GAPSO-optimized Fuzzy PID controller for electric-driven seeding. Through simulation experiments, speed accuracy experiments of seed metering device, and field experiments, it is verified that the proposed GAPSO is better than the traditional GA and PSO in the parameter tuning of the fuzzy controller; compared with the other two algorithms, the GAPSO fused by genetic screening mechanism has stronger global optimization ability under the same iteration time.; Macroscopically, the rotational speed accuracy of the seeding motor driven by the three different intelligent controllers of GAPSO-Fuzzy-PID is also higher, and the field experiments results are excellent, which provides the basic conditions for the high-density planting of crops.